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We study tunneling of charge carriers in single- and bilayer graphene. We propose an explanation 
for non-zero “magic angles” with 100% transmission for the case of symmetric potential barrier, 
as well as for their almost-survival for slightly asymmetric barrier in the bilayer graphene known 
previously from numerical simulations. Most importantly, we demonstrate that these magic angles 
are not protected in the case of bilayer and give an explicit example of a barrier with very small 
electron transmission probability for any angles. This means that one can lock charge carriers by 
a p-n-p (or n-p-n) junction without opening energy gap. This creates new opportunities for the 
construction of graphene transistors. 
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I. INTRODUCTION 

Klein tunneling, that is, transmission of massless Dirac 
fermions with a high probability throu^ potential bar¬ 
riers, whatever broad and high they arel^^l jg one of the 
key phenomena for graphene physics and technology. It 
protects conducting state of graphene with a high charge 
carrier mobility despite charge inhomogeneitiesSl at the 
same time, it does not allow to use the simplest con¬ 
struction of graphene transistor based on p-n-p (or n-p-n) 
junctions since such a device can never be lockecP. As 
a result, some tricky ways should be used, for example, 
tunneling transistor with vertical geometrji^ where elec¬ 
trons tunnel between two sheets of graphene. 

Full transmission for normally incident electron beam 
is symmetrically protected: since for massless Dirac 
electrons the propagation direction is intimately related 
to the direction of pseudospin, and the latter cannot 
be changed by an action of the electrostatic potential 
smooth at interatomic distances, the back scattering is 
completely forbidderfil before the discovery of graphene, 
this was noticed as an explanation of stability of con¬ 
ducting channels in carbon nanotubeJ^^. In the first 
calculation for the rectangular potential barrier in two- 
dimensional cas^, additional nonzero “magic angles” of 
the incidence were found with also full transmission; 
these magic angles were also found for the parabolic 
barriei® and associated to Fabry-Perot resonances known 
in optics. These additional resonan ces are not univer¬ 
sally protected. Semiclassical analysi^^E^^ has shown that 
these magic angles exist for symmetric potential barriers 
only whereas for generic one-dimensional barriers max¬ 
ima of transmission corresponding to the Fabry-Perot 
conditions are suppressed; moreover, this suppression is 


exponential in a formal semiclassical smallness parame¬ 
ter. 

Much less is known on the case of bilayer graphene 
where, in the simplest approximation, electron spectrum 
is massless but with parabolic touching instead of conical 
one and with nontri vial chir al properties of the charge 
carrier wave functiort- ^ l ^b2 | i3 | ^^30^ for the nor¬ 

mally incident electron beam the transmission probabil¬ 
ity is exponentially smalP^ this is a nice counterexam¬ 
ple to attempts to relate Klein tunneling in single layer 
graphene “just” with the gapless character of the en¬ 
ergy spectrum whereas the chiral properties of the wave 
functions are the most important. Existence of “magic 
angles” with full transmission in bilayer graphene has 
been demonstrated numerically for rectangular potential 
barrieil^ and for some smoothen shapes of the barrieJ^. 
In the last paper, it was claimed that for asymmetric po¬ 
tentials, contrary to the case of single layer, magic angles 
survive. As we will see this statement is not quite ac¬ 
curate. Chiral effects in penetration of charge carriers 
through potential barriers in bilayer graphene has been 
studied also in Refs. M and m However, the issue of 
stability of magic angles in the bilayer graphene is still 
unclear. They have no obvious symmetry protection, like 
100% transmission at zero incident angle for the single 
layer graphene. Apart from theoretical interest this is a 
question of potentially great practical importance: if it 
would be possible to create a potential barrier with small 
enough electron transmission probability for any angles 
this would open a way to build a conventional transistor 
based on p-n-p (n-p-n) junction in bilayer graphene with¬ 
out gap opening. In this paper we will show that this is, 
indeed, theoretically speaking, possible: the magic angles 
are not stable and electron transmission can be strongly 
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suppressed by a proper choice of the shape of the barrier. 


II. BASIC EQUATIONS AND THE 
FORMULATION OF THE PROBLEM 

A. Single-layer graphene 

Quantum mechanics of charge carr iers in graphene is 
governed by a massless Dirac equatiorP^ 

(I) 




Transmission probability Transmission probability 


(a) 


(b) 


H 


VpicTxPx + (^yPy) + U{x, y) 

( U{x,y) VF{Px-iPy)\ 
\VF{Px+iPy) u{x,y) J’ 


( 2 ) 


where E is an energy of the stationary state, ax-, ay are 
Pauli matrices, U{x,y) is the potential, Vf « c/300 is 
the Fermi velocity (c is the velocity of light), Px = —ihdx, 

Py = —ihdy. 

To study Klein tunneling, we restrict ourselves to the 
conventional case of a one-dimensional (i.e. given by a 
function depending on one variable U(x,y) = U{x)) po¬ 
tential barrier. A natural step is then to study the trans¬ 
mission and reflection probabilities as a function of the 
angle of incidence 9. As mentioned in the Introduction, 
the transmission probability at 6* = 0 is always 100%, ir¬ 
respective to the parameters of the potential. Numerical 
calculations show existence of 100% transmission at addi¬ 
tional magic angles for symmetric potential barriers and 
just maxima at some angles for the asymmetric one d^ l ^^ l 
To illustrate this, we show in Fig. our computational 
results obtained for the same potentials of n-p-n junc¬ 
tions which were considered in Refl3 that coincide up 
computer precision to those oP (Fig. 6). 

From a general point of view, the existence of the magic 
angles for symmetric potentials and their disappearance 
for generic ones is an unexpected property of the Dirac 
equation <§■ In Refsin and [TTI it was explained within a 
quite comp icated semiclassical theory. It is interesting 
to study its origin per se. Perfect transmission assumes 
that the reflection amplitude is zero; this means that a 
comp lex-valued function of one real argument, the angle 
of incidence, has nontrivial zeroes. We will give below a 
simple solution of this problem based only on symmetry 
properties of the Dirac Hamiltonian. 


B. Bilayer graphene 


The situation is different for the bilayer graphene. Fol¬ 
lowing ReflSj we will deal with the simplest effective 
Hamiltonian describing chiral particles with parabolic 
dispersion: 


( U{x,y) {p,,-ipyf/2m\ 
\{px + iPyY/‘2m U{x,y) J’ 
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FIG. 1. Numerical simulation of transmission probability for 
n-p-n junction in single layer graphene. Energy of a particle 
80 meV, height of the junction 200 meV. (a) Transmission 
probability for a symmetric potential of barrier width h = 
250 nm, for which n-p and p-n regions have characteristic 
widths Zi = Zs = 100 nm. (b) Transmission probability for an 
asymmetric potential with Zi = 150 nm and Z3 = 50 nm. (c) 
and (d): the corresponding potentials. Se^ (Fig. 6). 




where m « O.OSlme is the effective mass of the electron 
in the bilayer graphen^i^, We is the free electron mass. 
(This value is according to the latest experimental data; 
in early papers, the value m = 0.054me was used.) 

This Hamiltonian is not applied to the real bilayer 
graphene for very low energies {E < 1 0 meV ) where the 
effects of trigonal warping are essentiaP^E^ and for high 
enough energies {E > 200 meV) where the transition to 
four-band pictur^^^ is required. Anyway, Eq. ([^ is a 
new type of wave equation different from both nonrela- 
tivistic Schrodinger and Dirac equations and its study is 
by itself of significance for mathematical physics. 

Contrary to the case of single-layer, there is no more 
Klein tunneling at zero angle of incidence and no clear 
symmetry properties which would protect full transmis¬ 
sion at other angles, thus giving a hope to find a barrier 
that allows blocking of a current. 

Numerical experiments for an n-p-n junction in bi- 
layei!^ have shown the following: 

• nonzero magic angles are still present in a symmet¬ 
ric potential; 

• moreover (and most surprizingly!) they seemed to 
survive when one passes to a non-symmetric poten¬ 
tial, se^ (Fig. 7) and Fig. Hb) with our computa- 


( 3 ) 
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tional results. 

Both these effects thus were to be explained. 
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FIG. 2. Numerical simulations of a transmission through 
an n-p-n junction in bilayer graphene. Energy of a parti¬ 
cle 17 meV, height of the jnnction 50 meV. (a) Transmis¬ 
sion probability for a symmetric potential of barrier width 
h = 100 nm, with n-p and p-n regions of characteristic width 
li = h = 10 nm. (b) Transmission probability for an asym¬ 
metric potential with h — 20 nm and I 3 = 40 nm. (c) and 
(d): the corresponding potentials. 
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C. The specific goals 


In this paper we are answering all the above questions 
as well as obtaining some other results. Namely, we 

1. Obtain an equivalent condition for a magic angle for 
a symmetric potential. This condition is given by 
one real equation (with generically simple roots) . 
This is done for both single-layer (see Sec. Ill A) 


and bilayer (Sec. IIIB) graphene. In particular, 


the presence of magic angles for the symmetric n- 
p-n junctions (considered in Refs. [7l and fTTll cannot 
be removed by a small perturbation of the potential 
within the class of symmetric ones. 


2. Explain (see Sec. Ill) the seemingly magic angles 
for the nonsymmetric potentials (red line on Fig.j^. 
It turns out that the minimal reflection probabil¬ 
ity is very small (of order 10“^.. 10““^ for the po¬ 
tential studied in Ref.^ but is still nonzero (cf. 
Fig. [T]). The reflection/transmission probabilities 
for this potential turn out to be very close to the 
ones for a close symmetric potential (cf. Fig. |^. 


The latter possesses exact magic angles, that be¬ 
come a seemingly magic for a non-symmetric one. 


3. Find (see Sec. IV) an explicit example of a potential 
for which in some band of energies the transmission 
probability is less then 2 • 10“® for any angle of in¬ 
cidence. Interestingly enough, such a potential can 
be taken to be symmetric, in particular, showing 
that magic angles are not obliged to be present for 
an arbitrary symmetric potential U{x) for bilayer 
graphene. 


4. Provide (see Sec. VB) an approximation for the 


transmission probability p{0) around a peak, show¬ 
ing that such a peak has (approximately) a stan¬ 
dard Lorenz—Breit—Wigner form 


p{0) 


Po 


l + c{0- 0oy 


5. Provide (see Sec. VC) a “two-level” method of 


finding the transmission probabilities in bilayer 
graphene that removes the exponential growth 
problem. 


III. SYMMETRIC CASE: TRUE MAGIC 
ANGLES 


The key argument in both symmetric cases will be the 
following symmetry of the graphene equations: 


y^2{x,y)J yxif^(-x,-y)J- 


(4) 


It maps the solution of Q for a potential U{x,y) to 
the solution of Q for the potential turned by 180° 
U{—x,—y). Hence if U{x,y) = U{x), as we will assume 
through the rest of the paper, and U{—x) = U{x) (as we 
consider now the symmetric case), the map T sends the 
solutions of Q to the solutions of Q. Finally, it is easy 
to notice that T preserves the direction of a flat wave. 

It turns out that this symmetry of the problem reduces 
the number of independent equations. This is the same 
scenario that occurs, for insta nce, f or the equation de¬ 
scribing the Josephson junctioii ^^FS I 


A. Single-layer case 


For the sake of mathematical completeness, let us 
state the problem formally. As U = U{x) does not de¬ 
pend on y, the eigenfunctions can be tried in the form 
'i'{x,y) = 41 ( 0 ;, 0)e“^, where the corresponding eigen¬ 
value of the operator py is equal to ah. After denoting, 
by a slight abuse of notation, 'I'(a;) = ll'(a;, 0), the equa¬ 
tion Q becomes the ordinary differential equation 


f U{x) - E 
\-iVFh{dx - a) 


—lYph^dx + a)^ 
U{x) -E ) 


T'(a:) = 0, 


(5) 
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or, equivalently, 

/ ■ U(x)-E \ 

d^'S{x)=[ .ut)-E * 

In the domain U = 0 (that is, to the left or to the 
right of the barrier), the solution of ([^ is a linear combi¬ 
nation of left- and right-going waves iphix) = e~^^^VL 
and tpR{x) = respectively. Here ±.ik, where 

a? + k‘^ = , are the eigenvalues of the right hand- 

side operator in and 

( a + ik\ , f a — ik\ 

“ \iE/VFfi) “ [iE/VFh) 



are the corresponding eigenvectors. For a wave of energy 
E falling at the angle 9 we thus have 

E E 

a = —— sin 0, k = —— cos 9; 
ryf rYf 

we can also rewrite the eigenvectors as 


VR = 


E 

VrIt, 


and vl = 


E 

VrIt, 


,-i0 


( 7 ) 


To find the transmission and reflection probabilities for 
a given 9, one then looks for the solution of ([^ that is of 
the form 


^ ^ I '4’r{x) + r{9)tpR{x) on the left of the barrier, 
I t{9)'ipR{x) on the right of the barrier. 

( 8 ) 

Here r{9) and t{9) are respectively reflection and transi¬ 
tion amplitutes, and the corresponding probabilities are 
the squares of their absolute values. 

In a general quantum mechanical setting, these ampli¬ 
tudes are general complex numbers with \r{9)\^-\-\t{9)\^ = 
1; the “no-reflection” condition r = 0 on a complex num¬ 
ber r cannot be satisfied in a generic one-parametric fam¬ 
ily. It turns out, that in our particular situation of a 
symmetric potential, the numbers r{9) and t[9) satisfy 
an additional relation. Namely, we have the following 

Proposition 1. The function q{9) = where 

r{9) and t{9) are defined by Q, is real-valued for all 
values 9. 

Instead of the function q{9), it is more convenient to 
consider its differently normalized version: the function 
f(9) = , This real-valued function, on the one 

hand, satisfies \f{9)\ = |r(0)|, and hence its zeros are ex¬ 
actly the magic angles: see Fig. On the other hand, 
this real-valued function generically has simple zeroes 
(and as it takes values of both signs, its zeroes cannot 
be removed by a small perturbation). 


FIG. 3. Blue line: the transmission probability for the sym¬ 
metric potential from Fig. (a,c). Green line: real-valued 
analytic function f{6) such that |r(0)| = |/(^)|- Its zeroes 
correspond to magic angles. 


Proof. First note that the map T, given by Q, descends 
on the space of solutions of § , preserving a: 


T: 






(^^2(-x)\ 


( 9 ) 


Second, remark that the flow § preserves the bilinear 
antisymmetric form with the matrix 


Q = 



that is, for any two solutions the value 

Q(' 0 i,' 02 ) = (“t” means transposed) does 

not depend on x (and thus Q becomes a well-defined bi¬ 
linear form on the space of solutions). The easiest way 
to see it is to say that the matrix of the flow ([^ has zero 
trace, and hence the flow map from a: = xi to a; = 0:2 has 
determinant one, while ' 0 iQ '02 = det{ipi 'ip 2 )- 

It is also immediate to check that Q{T'tpi,T'tp 2 ) = 
Q(i/’ 2 , f/’i), that can be re-stated as 


Q{fii,Tfi2) = Q{tf2,T'ilJi). 

This remark easily implies the following observation: for 
any solution ip of the value Q{'tp,T'ip) is purely real. 

Now, let us apply this remark to the solution 41 = 
j^d/R, where 4'^ is of the form Q. On one hand, 
(5(41, rtk) is a real number. On the other hand, cal¬ 
culating it at any point a;o on the left of the barrier, we 
get 


Q(4',T4') = Q{j^fi’Rixo)-\- 


r{0) 

t{9) 


Q{vl,Tvr) = 


rjS) 

m 

rjd) 

m 


'(Pl{xo), (rV’fl)(xo)) = 


• (e 


- 2 ie 


1 ) 


E^ 

w 






















where 
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where we have used Q for the last equality. Then, 

2aE 
Vph' 




rp2 rp2 

= —• 2i sin0 • 






= —le 


and hence 


Q(vE-,TvI/) =-g(0) 


2aE 

V^' 


The left hand side of the last equality is real, and thus 
so is q{9). 

□ 

To conclude this paragraph, let us restate the condi¬ 
tion for an angle being magic, r{0) = 0, in two different 
ways. First, due to the proof of Prop, [^it is equivalent 
to = 0; calculation of this form at a; = 0 

gives us 

Q{^r,TM>r) = |^fl,i(0)p - |4'fl.2(0)p. 

Second, the form Q is non-degenerate, so Q('*/'i, V' 2 ) = 0 
if and only if ijji and '!/'2 are proportional. We finally get 

Proposition 2. The angle 9 is magic if and only if 

l'^'fl,i(o)| = I'i'i?, 2 ( 0 ) 1 , 

and if and only ifT'^R is proportional to 

B. Bilayer case 

In the same way as in the single-layer case, we are 
considering the solutions of the form 'I'(x, y) = 'I'(a;)e“^. 
Then equation 0 becomes 


/ U{x) E = 0 

[-tid.-a? U{x)-E ) 


( 10 ) 


or, introducing new variables 

Uix) = ^(Uix) - E), ^i{x) = - a)^i{x), 

^2ix) = (d^ + a)^2ix), ^ 


we can reduce (101 to the following equation: 

( 




a 1 0 0 \ 

0 a U(x) 0 

0 0 -a 1 

\U{x) 0 0 -aj 




( 12 ) 


on a complex 4-vector 

^(cc) = (^'i(a;),$i(a:),4'2(a;),$2(a;))*. 

In the same way as for the single-layer case, in the 
domain {7 = 0, equation (12) has the solutions of the 


form (for E > 0, otherwise some signs will be different) 

tpR{x) = e^^'^VR, tpL{x) = 


/ a + ik '\ 


VR = 


2m TT 
h? ^ 

— {a — ik) 

) 


/ a — ik \ 


and vl = 


— {a + ik) 

) 


(13) 


and E k^ = '^E- The angle of incidence 9 is now 
related to o by a = ^^E sin 9 and k = ^^E cos 9. 
Though, the equation (12) has also in the domain U = 


0 the solutions of the form 


where 


tf+{x) = e^^v+, 

d-{x) 

1 

< 

1 

II 


a + A ^ 

2m 77’ 
h? ^ 

a — X 


2m rp 

Q, X 



, V- = 


/ 9 2 to 

-¥e) 


[-¥e) 


(14) 


v+ = 


To find the transmission and reflection probabilities 
one now has to consider the solution of the form 


^'fl(a;) = 


J dnix) + r(9)ipL(x) + aiip+ix), 
^t(9)dR(x) + a2d-(x), 


where suppC/ C [—a:o,a;o]. 

We then have the following 


X < —xo 
X > Xo 

(15) 


Proposition 3. The function q{9) = where r{9) 


andt{9) are defined by (151, is real-valued for all values 9. 


Once again, instead of the function q(9) we can con- 

which is more con- 




sider the function f(9) = , , , 

^ ’ Vl+9"(S) 

venient due to the relation |r(0)| = \f{9)\. Zeroes of 
f{9) correspond to the magic angles, and are generically 
simple (see Fig. |^. 

Proof. In the same way as in the single-layer case, the 
map T descends on the space of solutions of (12), pre¬ 
serving a, and becomes: 


T: 


/^'i(a;)N 

d>i(a;) 

^' 2 ( 2 :) 

Vd'2(a;)y 




/ ^^(-a;) \ 
^2(-a;) 

^j_{-x) 

\-^i{-x)J 


(16) 


Also in the same way as before, we note that the flow 
of the equation (12) preserves a bilinear antisymmetric 


form. Namely, note first that is preserves a sesquilinear 
form with the matrix 


Q = 


/ 0 0 0 IN 

0 0-10 
0 10 0 
\-i 0 0 Oy 
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FIG. 4. Blue line: the transmission probability for the 
symmetric potential from Fig. [^(a,c). Green line: real-valued 
analytic function f{ 6 ) such that |r(0)| = \f{S)\- Again, its 
zeroes correspond to magic angles. 


and that due to the invariance of Q by the flow the Q- 
product Q{vi,V 2 ) on two eigenvectors vi,V 2 can be non¬ 
zero only if the sum of corresponding eigenvalues van¬ 
ishes. □ 


Take now 

^ f X < - 3 ^ 0 , 

X>Xo. 

Then, Q(1E', T{'^)) is a real number. Note that at a: < xq, 
we have (T'I')(a::) = —ipfi{x) + {a 2 /t{ 9 ))ip+{x)), and hence 


Q('F,rw) 


rjd) 

m 


Q{vl,Tvr) = -Aik ■ 


2mE r{9) 

~W~ ' 


(The expression in the middle is the only term that does 
not vanish.) As the expression in the left hand side is 
real, so is q{9) = □ 


This can be checked by explicit computation; though, 
there is a simple physical interpretation for it. Namely, 
i$*(a;i)Q<i)(a:i) measures the current density through 
the section x = xi (see, e.g.^ (Eqs. (136)-(137))), so 
the conservation of the form is merely the conservation 
of number of particles. 

At the same time, both matrices of the flow and of the 
form Q are purely real. Hence, the flow also preserves 
an anisymmetric bilinear form with the same matrix Q. 
Hence (5(-, •) is a well-defined antisymmetric form on the 


space of solutions of ( 12 |: the value ('k(i)(a:))*( 5 lk( 2 )(x) 
does not depend on the choice of the point x. 

Finally, it is easy to see that we have the following 


Lemma 1. For any two solutions of (12) one has 


Q{i’i,Tilj2) = Q(V'2,F^/>i). 


In particular, Q{ipi, Tipi) is always real. 

For the eigenvectors of the matrix of the flow { 12) 
at zero potential, the only pairs giving nonzero product 
Q{ipi,T^ 2 ) are 


IV. TRANSMISSION-BLOCKING EXAMPLE 

The above arguments show that in the symmetric case, 
if magic angles were present for some potential U, they 
usually can not be removed by its small perturbation. 
Though, these arguments do not imply that the magic 
angles should be present for any symmetric potential: 
indeed, a real function q{9) generically is not obliged to 
have real roots. 

Indeed, one can construct an example of a symmetric 
barrier, the probability of transmission through which 
is quite small. The potential that we construct is 
a sufficiently quickly oscillating one, being a series of 
four n-p-n barriers. Namely, fix the height of barri¬ 
ers Uq = 50meH and the n-p and p-n junction widths 
h = h = ■ ‘' = hs = lOnm and the pairwise distances 
and widths I 2 = I 4 = ■ ■ ■ = li 4 = lOnm. Then, take 

So = 0, Xi+i=Xi + li, i = l,...,15 
and define 


Q{vr,Tvl) = -Q{vl,Tvr) = Aik ■ 

and 

Q{v+,Tv+) = -Q{v_,Tv_) = -AX ■ 

Proof. The first part is immediate; the equalities in the 
second part can be checked by an explicit computation. 
Finally, to check that all the other likewise products van¬ 
ish, note, that 

T{vr) = -VR, T{vl) = -VL, T{v+) = V-, T{v-) = v+, 


' 0 , if s < So or s > Si 5 or s £ [s 4 i_i,S 4 i], 

Uq, if S £ [s42-|_i , S4j_|-2], 

C/o • 5 (1 + tanh( 10 (s - , 

' if S £ [s4i,S4j+l], 

Uq- ^ (l- tanh(10(s - , 

, if S £ [s4i+2, S4i+3] 

(17) 

(see Fig. compare withP (Eq. (135))). 

The fast oscillations of the potential U (s) prevent the 
appearances of “resonances” between the junctions, and 
allow to block the transmission in a band of energies 
sufficiently close to C/o/2. Namely, in the energy band 
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FIG. 5. Fast-oscillating potential U{x). 


from 20 meV to 30 meV, the transmission probability 
p{9) = \t{9)\^ for any angle of incidence 9 does not ex¬ 
ceed 2 • 10“®: see Fig. 

Further increase in number of barriers allows further 
reducing of transmission probability in this band (or of 
slight widening of the band where a given upper estimate 
on |t(0)P holds). We will discuss in the next section a 
reason of why a blocking potential should be looked for 
among fast-oscillating ones. We conclude this section 
by noticing that other fast-oscillating potentials exhibit 
similar behavior. For instance, considering the potential 

\{l-cos{'Kx/2li)), a:G[0,16Zi] 

10 otherwise, 

we get a similar behavior for the maximal transmission 
probability: see Fig. top. 


V. BILAYER GRAPHENE: GENERIC 
ONE-DIMENSIONAL POTENTIAL 

A. Approximate magic angles 

We start by considering the example that was studied 
in ReflTl an asymmetric n-p-n junction. First, note that 
the corresponding angles are not exactly magic: though 
they seem to be such, precise computations show that the 
local minimal values of |r(0)| in these cases are « 0.029 
and « 0.053 at the angles 9 « 30° and 9 « 65° re¬ 
spectively (and hence the local maxima of the trans¬ 
mission probability |t(0)F are respectively « 0.9991 and 
« 0.997); see Fig.[^ 

Thus the correct mathematical question is not to ex¬ 
plain the exact equality but to explain why the minimal 
value of |r( 0 )| is so small. 

The first idea here would be to compare this asymmet¬ 
ric potential to a close symmetric one. Namely, in addi- 


(a) 
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(b) 

FIG. 6 . (a) Green line: maximal transmission proba¬ 

bility Pmax = max 0 p{0) for the potential (171 as a func¬ 
tion of the given energy E of the wave. Note that in the 
band E=20..30 meV the transmission probability does not 
exceed 2 - 10“®. Blue line: maximal transmission probabil¬ 
ity for the cosine-oscillating potential ( |18[ ). Note, that in the 
same energy band the transmission probability is still very 
small (it does not ex ceed 10 “®). (b) transmission probability 
for the potential (171 as a function of the angle 6 for energies 


E = 20, 22.5, 25, 27.5, 30 meV. 


tion to the asymmetric potential with junction widths 
= 20 nm and (3 = 40 nm and flat p-part of width 
I 2 = 100 nm, consider a symmetric one with the junc¬ 
tion widths /i = (3 = 20 nm and flat p-part of width 
I 2 = 110 nm (so that the middle of the p-n junction does 
not move). It turns out that the transition probabili¬ 
ties for these barriers are quite close to each other (see 
Fig. § red and green lines). At the same time, latter 
potential possesses exact magic angles due to the same 
arguments as in Sec. |IIIB[ Given this, one could argue 
that the reason for almost-magic angles is just that the 
solutions of two equations, the one for the symmetric po- 



































Inc. angle (degree) 

(a) 



Inc. angle (degree) 

(b) 


FIG. 7. Zoom in: absolute value of the reflection amplitude 
|r(6')| for n-p-n junction in bilayer graphene near magic angles. 
Settings are the same as in Figure]^ (b,d). 



Incidence angle (degree) 


tential and for the asymmetric one, are sufficiently close 
to each other. 

However, a closer examination of the graph on Fig. 
shows that this argument alone cannot be a satisfac¬ 
tory explanation. Indeed, one notices that the peaks are 
shifted with respect to each other, so that the value of 
one of the functions at the maximum of the other one is 
quite far from 1, much farther than its maximal value. 
Hence, additional arguments are required for an expla¬ 
nation here. 


In fact, note that altering both junction widths, that 
is, considering the symmetric potential with h = I 3 = 
30 nm, I 2 = 100 nm, one gets a much better approxi¬ 
mation for the transmission probability: see Fig. red 
and blue lines. Though, a priori it is not clear why such 
an approximation (contrary to the one-side modification) 
gives so precisely the approximately magic angles. We 
will explain it at the end of this subsection. 


To provide a complete explanation for the almost- 
magic angles effect, we will approximate the problem of 
crossing of an n-p-n barrier as a sequence of two inde¬ 
pendent crossings, of an n-p and of a p-n barriers respec¬ 
tively. The error in such an approximation will be almost 
neglectable, thus reducing our question to the study of in¬ 
dividual crossings. Finally, an additional effect, appear¬ 
ing in the bilayer graphene (contrary to the single-layer) 
is in the core of the explanation here. 


Namely, under the width I 2 flat part of the barrier, the 
equation (12) again has constant coefficients, and can be 
interpreted as zero-potential equation for a wave of en¬ 
ergy £)[/ :=E — Uo. Provided that \Eu\ > |a| = |i?sin0|, 
the wave-type solutions of (12) in this domain can be 
written in the form 




where kjj = \JE^ — a?. Here and v'j^ are the corre¬ 
sponding eigenvectors of the matrix of the system, given 
by (131 for the energy Eu, up to some sign changes due 


FIG. 8. Transmission probability for n-p-n junction in bi¬ 
layer graphene for asymmetric 20-100-40 nm (red lines), sym¬ 
metric 20-110-20 nm (green lines) and 30-100-30 nm (blue 
lines) barriers. 


to the inequality Eu < 0: 



/ a + iku\ 


(a - iku\ 

I’R = 

2m 77' 
h? 

a — iku 

II 

2m u 

a -I- iku 


\-¥Eu) 


K-^Eu) 


Also, as earlier, one also finds in this domain solutions of 
the form 


where Xu = ^/E^^+lE, and v± are the associated eigen¬ 
vectors (again with a sign change with respect to ([T4|)): 


/ a + Xu \ 


/ a - Xu \ 

2m u 

n? Eu 

, = 

2m 77' 
h? 

-{a - Xu) 

— (a + Xu) 

V -^Eu ; 


[ -¥Eu y 


Now, related to the problem of describing an n-p-n 
barrier, one can consider the problem of describing an 
n-p transmission, given by a potential Un-p(x) that is 
identically 0 on the left of the barrier and identically Uq 
on the right of it. For a barrier of characteristic width I, 
analogously to the n-p-n barrier, we can take 


U{x) 


0, if cc < 0, 

Uo, ii X > I, 

(7o • I (1 -I- tanh(10(a: - ‘^)/l)) , 
if a; G [0, 1], 


( 21 ) 


(compare witfP (Eq. (135))). A physically meaningful 
solution then should be of the form 


axiiR{x) E a2'4>L{x) E a3^/>+(x), 
to the left of the barrier, 
ai4>R{x) + a^ipuix) + aG(j)-{x), 
to the right of the barrier. 


( 22 ) 
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The coefficients 01,02 of the wave component to the 
left of the barrier and the coefficients 04 , 05 to the right 
of it are then related by a transmission matrix A„_p = 


where Ra stays for the rotation at angle a. 

Next, consider the p-n transmission problem. To do so, 
note, that the application of T sends it to the solution of 
the form 


Oi 

02 



Note now, that for any coordinates on the the 2- 
dimensional space of physical solutions, the restriction 
of the bilinear antisymmetric form Q on this plane is 
proportional to the determinant (area) form in these co¬ 
ordinates. Hence, considering two different systems of 
coordinates ( 01 , 02 ) and ( 04 , 05 ), we see that the coeffi¬ 
cient of proportionality is equal to Q{tj}R, '>Pl) = 4ffc- 
and Q{4>r, (jih) = iiku ■ ^Ejj respectively, and thus the 
determinant of the matrix An-p, relating these coordi¬ 
nates, is equal to 


det Aji_p 


Aiku ■ ^Eu 
^ik ■ ^E 


kuEjj 

kE 


Now, pass from the solutions ipRjipL in [/ = 0 domain 
and from in the U = Uq one to the properly nor¬ 

malized “sine-cosine” solutions 


'Pcos 


VkE 

2 


{ipR+IpL), 


'Psin • — 


2 i 


{ipR - tpL) 


^ p—n (^) 


aiipR{x) + a2'ipL{x) + O3'0_(a;), 
to the right of the barrier, 
a44>R{x) + a^cpLix) + afi4>+{x), 
to the left of the barrier 


(26) 


for the new potential Up-n(x) = Un-pi—x) (note that 
this potential will be supported on [—Z,0]). For the cor¬ 
responding transmission matrix, we then have 

Ap-n = A^_p. (27) 

In the same way as before, we pass to the sine-cosine 
bases, thus obtaining a real determinant 1 matrix 


A — 

-^p—n — 


yJ—kjjEu 

VkE 


(-i i) (i 


Joining (|23|), (|27| and (|28|), and using the reality of ma- 

(29) 


trices An-p, Ap-n, we get 


A — (1 0 \ 4-1 0 'i 

^p-n — f 0 -1 j ^n-p f 0 -1 j ■ 


in the domain U = 0 and to 
\U-kuEu 


ipR-Eph), Psin ■ — 


\/—kuEu 

2 i 


{(pR-pL) 


in the domain U = Uq. The advantage of these solutions 
is that as they are purely real, the transmission matrix 
in these coordinates 


4 

-^n—p 


'/kE / 1 
/—kuEjj V 21 




(23) 


is also purely real, and due to the choice of the normal¬ 
ization it is of determinant 1 . ^ 

Recall now that a real area-preserving matrix A admits 
singular value decomposition: it can be represented as a 
product of a rotation matrix, a diagonal matrix ^ ^ 1 j , 

and another rotation matrix. The reflection and trans¬ 
mission coefficients, associated to such a matrix then sat¬ 
isfy 


\t\ 


2 

i — 1 ’ 
+ 




(24) 


In particular, for an angle to be (approximately) magic, 
the corresponding real matrix should be an (approxi¬ 
mate) rotation. 

Write for a transmission problem through an n-p bar¬ 
rier of characteristic width I and for the angle of incidence 

e 


An—p — An—pil, ^) — Ra{l,G) 


(Ai,e) 0 

V 0 


^ 0 ( 1 , 6 ), 

(25) 


The singular decomposition ([2^ thus gives us 


4 — (i 0 \ Tf (Uifi) 0 ^ Y 

— (0-1 j0 ^ 

X R-a{l,e) ( 0 -°1 ) = Rp(l,e) ( ^^' 0 ^ ) Ra{l,0) ■ 

(30) 


Finally, let us return back to an n-p-n junction. Note, 
that if the width I 2 of the flat p-part of the barrier is 
sufficiently large, the transmission through the barrier 
can be approximated as a sequence of an n-p and p-n 
junctions (of widths R and I 3 respectively). Indeed, the 
corresponding physical solutions are of the form 


'kn—p—n(^) — ^ 


aitpRix) + a2ipL{x) + aztp+ix), 
to the left of the barrier, 

UipRix) + a^p/x) + aep-ix) + ayp+ix), 
in the p-zone, 

asi>Rix) + a^ipLix) + aioZ-ix), 
to the right of the barrier. 


The component aTp+ix) is of order at most 1 near the p-n 
part, and hence is of order near the n-p transition, 

which is neglectable (it is less than 10 “^ for the barrier 
and angles described on Figs.[^|^. The same applies to 
aQp-ix). Thus, with very high accuracy one has 
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where the transmission matrix corresponds to the 

p-n barrier, shifted by T = /i + Z 2 + ^3 to the right from 
the position [—^ 3 , 0 ] at which it was studied earlier. For 
the full transition matrix we thus have 


A 


n—p—n 


Aji—pA 


(!') 

p—n' 


Hence, the same holds for the transmission matrix in 
the sine-cosine basis: 


A -A A^‘'^ 

-^n—p—n ~ -^n—p-^p—n' 


(31) 


Taking into account that 


what is in perfect agreement with |r(0i)| « 0.0299. For 
02 ~ 65.34° one has 

log/ri( 6 » 2 ) « 1.5815, logAt 2 ( 6 ' 2 ) « 1.6339, 

S( 02 ) = I log/ri( 6 » 2 ) - log^ 2 ( 02)1 « 0.0524, 

what is also in perfect agreement with |r( 02 )| ~ 0.0524. 

Finally, the above description also explains why the 30- 
100-30 nm potential was such a good approximation for 
the 20-100-40 nm one (see Fig. [^. Indeed, the equation 
for the magic angle is based on the angles /3(0, li),j3{9, 13 ), 
and we have with quite high precision 


^p-n ~ R-kjjl' Ap-nR, 


kl' 


and substituting into (31) the singular decomposi¬ 
tions (25), (30), we obtain for the total n-p-n transmis¬ 
sion matrix 


An—p—n{0') ~ R* 


Qi(e) 


where 


/ Ml 

V ' 


X R 


( 9 ) 0 

0 

-1 


Rff 


{ 9 )- 


■02 {0) — kul' 


0 


V 0 p.2{0) 


R. 


'Oi 2 {&) — kl' 5 


(32) 


0 i{S) = 

ai{ 0 ) = a{ 0 ,li), 

f3ii9) = I3i0,h), 


9-2(0) = iJ.{ 9 ,l:i), 
a2(9) = a(9,h), 
P2(9) = P(9M). 


This matrix is natural to expect to be closest to a ro¬ 
tation for angles 9 when the composition of rotations in 
the middle, Rp^(^ 0 )+ 02 ( 9 )-kui'^ is a rotation by an integer 


multiple of tt (compare with (33) and (34) for in Sec. VB 
below). For such angle 9o, one has 


+ 13(11, k) 


mo. 


To conclude, we note that the mechanism described 
in this section is exactly the one that one would like to 
avoid while looking for a transmission-blocking potential. 
Hence, it seems quite natural to avoid long flat parts 
(as they are likely to “cancel”, for some angles, what 
happens before and after them), and hence consider a 
fast-oscillating potential — as it was done in Sec. IIV 


B. Peaks for n-p-n barriers 


The description above can be used to describe the 
shape of a peak of the transmission probability. Namely, 
for a matrix A G SL{2,M.), the corresponding inertia 
coefficients ^ and can be found from the relation 
trAA* = thus implying that 


p(9) = m\^ 


4 

2 -f tr AA* ’ 


(33) 


An-p-niO) « ±R, 


ni9) 


p:[A9)p2{9) 


R. 


■a 2 id) — kl' • 


Now, a final (and key) remark is that for the bilayer 
graphene and relatively short n-p barriers, the value 
p,{9,1) does depend on the angle of incidence, but de¬ 
pends very slightly on the width of the barrier, while the 
latter stays sufficiently small: see Fig. |9| A plausible 
explanation for this will be discussed in Sec. |V C| 

Due to the this observation, the above matrix is indeed 
very close to the rotation one, as fi{9, li) « fj,(9, 13 ). More 
precisely, for such 0 q, substituting fi = 9i(9[j)/^ 2 ( 00 ) 
into (24), one has 


where A = An-p-n is the corresponding transmission 


matrix. From the approximation (32), one gets 




' pI(9) 0 

0 p-\9) 


^ R—j{ 9 ) 1 


where 7 ( 0 ) = (3i{9) + 132 ( 0 ) —kul'■ An explicit calculation 
then gives 


trAA* = {fj.l{9)pl(9) + ( 9 )^ 2 '^ (9)) cos'^ j(9)+ 

+ (91(0)92^(0) + (9)nl(9)) sin^ -f(9) = 

= ci(9)+c2(9)cos‘^j(9), (34) 


k(^o)| = ^ J = tanh(5(6>o), 

9 + 9 

where (5(0) = | log/ri(0) —log^2(0)|. This scenario indeed 
holds for both approximate magic angles illustrated on 
Fig-d right. For instance, for an approximate magic 
angle 0i « 30.65° we have 

log/ii(0i) « 1.1641, log/X 2 ( 0 i) « 1.1343, 

S(9i) = I logpi(9i) - log^2(0i)| « 0.0298, 


where 


c,(9)=pI(9)p^\9)+p^\9)pI(9), 
C2(9) = (pI(9)-p^\9)){pI(9)-p^^{9)). 


A peak of the transmission probability corresponds to 
a local minimum of the denominator of (331. Due to the 
wide Z 2 -part, the angle 


j(0) = Pi(9) + P2(9) - kul' 
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0° 10° 20° 30° 40° 50° 60° 70° 80° 


Incidence angle (degree) 



Incidence angle (degree) 
(b) 


FIG. 9. The function log ^{9,1) for a wave of energy E = 
17 meV, crossing on bilayer graphene an n-p barrier of height 
V — 50 meV (same energy and height as for Fig. [^. (a) 
the dependence on the width I for (almost) magic angles 6 ~ 
30.64° (blue line), 9 rs 65.34° (red) and for 9 — 45° (green), 
(b) the dependence on the angle of incidence 9, for widths I = 
20 nm (blue line), 30 nm (green), 40 nm (red). 


changes much faster, than do, so such a peak Oq 

(almost) corresponds to 7 (^ 0 ) ~ Approximating 

cos^( 7 (d)) « c{0 — 0o)^, where c = (^ 7 ^)^, one gets for 
the shape of the peak the Lorenz-Breit-Wigner form: 


p{0) 


4 

2 + tr AA* 


P%) 

l + c'(0-0o)2' 


Such an approximation turns out to be quite precise: see 
Fig. for such an approximation for the peaks corre¬ 
sponding to the n-p-n barrier discussed irP (as well as in 
Sec.lHBl). 


FIG. 10. Blue (solid) line: exact transmission probability 
for the symmetric n-p-n junction. Red (dashed bold) line: 
Lorenz-Breit-Wigner approximation p{9) « • 


C. Computational algorithms 


We conclude this paragraph with specifying a compu¬ 
tation method for hnding the transmission probabilities 
for long barriers. Namely, a straightforward method of 
computation includes numerically solving the differential 
equation (12 1 through the barrier. Then, finding a linear 


combination 41^ of solutions that starts on the right of 
the barrier with ipR and with ip- that would have no ex¬ 
ponentially growing component on the left of the barrier. 
Though, if the barrier is sufficiently long, the solution 
starting with a given initial value has a tendency to grow 
exponentially with the width of the barrier. 


Instead, consider the solution 4>_(a;) to (12 1 that coin¬ 


cides with ip- on the right of the barrier. Note, that we 
do not need to know the solution itself, but only up to 
the proportionality: its only role will be to be added to a 
linear combination to remove the exponential growth on 
the left of the barrier. Hence, instead of finding 4>_(a;) 
itself (which is most natural to expect to grow exponen¬ 
tially when one passes to the left of the barrier), we can 
look for a unit vector-valued function (x) that is pro¬ 
portional to it at any point x, given by the normalization 


4*° (a;) = 4'_(a:)/|'I'_(x)|. 


(35) 


The latter can be found either by normalizing the solu¬ 
tion on each step, or by solving a differential equation for 
it, 


a,,4-° = ^4-° - (A4'‘).,4'° )4')!., 


(36) 


where A be the matrix of ( [T^ and (•, •) is the usual scalar 
product. 

Second, to find the reflection and transmission coeffi¬ 
cients, we are looking for a solution 'Ffl; to (121 that is 
of the form (15). Again, it is natural to expect that the 
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solution, starting with ipn on the right of the barrier, will 
grow exponentially. 

But the solution we are looking for is anyway a 
combination of the solution starting with ipn and the 
one starting with ip-. So while solving the equation (12) 
right-to-left, starting with ijjji outside the barrier, we can 
at any point x safely add (a;) with any coefficient: this 
keeps us in the same space of linear combinations. Thus, 
we consider the component of ipR that is orthogonal to 
'k'i, that is, 


ip{x) := iPr{x) - {ipR{x),'if°_{x))'S°_{x). (37) 


It obeys the following modification of the equa¬ 
tion (fl^: 


d^iP = AiP-]^{{A + A*)ip, vI/0_)d/0 . 


Such an orthogonalization removes the growth associ¬ 
ated to the highest (non-physical) eigenvalue of A, leav¬ 
ing only the expansion and contraction associated to the 
physical solutions themselves. 

Finally, after arriving to the point (—xq) on left side of 
the barrier, one further modifies the solution ip on the left 
of the barrier by adding again tk)). with such a coefficient 


Q{ip{-xo),v+) 

Q{^°_{-xo),v+) 


that the obtained combination ip := ip + /3'k^ does not 
have a component along ip-. Decomposing ip{—xo) in 
the base of eigenvectors vr,vl,v+,V-, we find and 

as coefficients before vr and vr respectively. 

The above algorithm (or its slight modification) seem 
also to be an explanation of the almost-constancy of the 
inertia coefficients /i(6*), that are (as it was discussed in 
Sec. V A) in the background of the appearance of almost- 
magic angles in the bilayer case. 

Note first, that we can construct the space of physically 
interesting solutions of (12) in the following way. Any 
bounded (“physical”) solution tk of (12) is Q-orthogonal 
to both solutions and 'k- that coincide respectively 
with '(/)_|_ to the left of the barrier and with ip- to the 
right of it. Indeed, Q-products (5('k,tk+) and (5('k,4'_) 
give (after normalization by Q{ip-,ip+)) respectively the 
coefficients before ip- in the decomposition to the left of 
the barrier and before 1 p.^. in the decomposition to the 
right of it. On the other hand, for a generic 6, E, U the 
solutions 'k_|_ and 'k_ are linearly independent, and hence 
a Q-orthogonal complement to the 2-plane {a4'+ -1-/3^'-} 
is exactly the 2-plane of physical solutions. 

Now, instead of considering functions 'k- and d>_|_, we 
can consider their normalized versions 


VkO (x) = d/+(x)/|^+(x)|, ^^-{x) = 'k_(x)/|vk_(x)|. 

As before, they can be constructed by solving a differ¬ 
ential equation left-to-right and right-to-left respectively. 


Note that it is natural to expect such a construction to be 
stable under small perturbations: a long-time map con¬ 
tracts almost all the directions in a small neighborhood 
of the most expanded image. 

Once these two functions are found, at any point 
X we know the 2-plane of physical solutions as the 
Q-orthogonal complement to the 2-plane {a'k(J_(x) -I- 
(x)}. 

Finally, knowing such plane, generically, we can recon¬ 
struct the flow without using the potential U explicitly. 
Indeed, the derivatives of the first and of the third coordi¬ 
nate in (121 do not include U, and at the same time, the 
first and the third coordinate generically provide a sys¬ 
tem of coordinates on the 2-plane of physical solutions. 
Thus, the system of differential equations on the first and 
third coordinate can be written only using tk^. 

The same technique applies also to the study of n-p 
transitions, with the only difference that on the right of 
the barrier we are taking the function instead of ^ 
that is obtained as <!>*(. = <I>_/|<k_|, where <I>_ coincides 
with <p- on the right of the barrier. 

At the same time, it is natural to expect that 'k)), and 
will not change too much under small changes of the 
potential U defining the n-p transition. The reason for 
that is that these functions are solutions to 1st order dif¬ 
ferential equations, defined by U (so even a discontinuity 
of U will result in a smooth behavior of and 
At the same time, in a zone where two different poten¬ 
tials U coincide, these functions (roughly speaking, cor¬ 
responding to the fastest-growing direction) quickly be¬ 
come aligned, so the effect of any perturbation should be 
sufficiently local. 

This shows the difference between the single- and bi¬ 
layer graphene. For a single-layer graphene, we have a 
two-dimensional differential equation directly involving 
the potential U. While for the bilayer graphene, we have 
an additional “integration”: the solutions associated 
to these potentials, are very close to each other (and the 
same holds for the solutions $° ); they mostly coincide, 
and differ only a little between the two potentials. And 
as it is these functions that define a new 2-dimensional 
differential equation, the difference between the trans¬ 
mission matrices An-p comes from the “integration” of 
difference between them, and should be even smaller. 


VI. CONCLUSIONS 

To summarize, we have considered chiral tunneling 
through one-dimensional potential barriers for the cases 
of both single-layer and bilayer graphene. We have 
proven that in both cases for symmetric barriers magic 
angles with 100% transmission can be found from one 
real equation. For the case of bilayer, this equation does 
not necessarily have real solutions and we have presented 
examples of the potential barrier with a very low trans¬ 
mission at any angles in a restricted energy range. This 
opens a way to build the conventional p-n-p (or n-p-n) 
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transistor from bilayer graphene without opening the en¬ 
ergy gap. This is also important conceptually since it 
gives a clear counterexample to an opinion that the Klein 
tunneling in single-layer graphene is due to gapless char¬ 
acter of the spectrum; actually, it is due to a special chiral 
character of electronic states there. We have also pre¬ 
sented arguments explaining why for the case of bilayer 
graphene the difference between symmetric and antisym¬ 
metric barriers as not as dramatic as for the single-layer 
one: whereas for the latter case asymmetric shape of the 
barrier results in an exponential suppression of transmis¬ 
sion at non-zero magic angles, in bilayer graphene they 
are robust, in a sense that the transmission probability 


remains close to 100%. 
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